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Magnetic properties of Fe, Co, and Ni at finite temperatures have been investigated on the basis of 
the first-principles dynamical CPA (Coherent Potential Approximation) combined with the LDA (Local 
Density Approximation) + U Hamiltonian in the Tight-Binding Linear Muffintin Orbital (TB-LMTO) 
representation. The Hamiltonian includes the transverse spin fluctuation terms. Numerical calculations 
have been performed within the harmonic approximation with 4th-order dynamical corrections. Calculated 
single-particle densities of states in the ferromagnetic state indicate that the dynamical efi'ects reduce the 
exchange splitting, suppress the band width of the quasi-particle state, and causes incoherent excitations 
corresponding the 6 eV satellites. Results of the magnetization vs temperature curves, paramagnetic spin 
susceptibilities, and the amplitudes of local moments are presented. Calculated Curie temperatures (Tc) are 
reported to be 1930K for Fe, 2550K for Co, and 620K for Ni; Tc for Fe and Co are overestimated by a factor 
of 1.8, while Tc in Ni agrees with the experimental result. Effective Bohr magneton numbers calculated 
from the inverse susceptibilities are 3.0 /ib (Fe), 3.0 /ib (Co), and 1.6 /ib (Ni), being in agreement with 
the experimental ones. Overestimate of Tc in Fe and Co is attributed to the neglects of the higher-order 
dynamical effects as well as the magnetic short range order. 

KEYWORDS: dynamical CPA, metallic magnetism, iron, cobalt, nickel. Curie temperature, effective Bohr mag- 
neton number, excitation spectra 
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1. Introduction 

The 3d transition metals, Fe, Co, and Ni are well- 
known to show the ferromagnetism with rather high 
Curie temperatures 1040 K, 1388 K, and 630 K.^^^) 
Magnetic properties of these materials are character- 
ized by the itinerant as well as local-moment behav- 
iors.''-' According to the angle resolved photoemission 
spectroscopy,^'^-' quasiparticle d bands with the Fermi 
surface are observed in these transition metals, and the 
noninteger ground-state magnetizations per atom in unit 
of the Bohr magneton number (/ie) are found;^'*-' 2.2 
HB (Fe), 1.72 (Co), and 0.62 (Ni). These re- 
sults indicate that the Sd electrons are itinerant. On 
the other hand, the Curie- Weiss law in the susceptibil- 
ity with atomic effective Bohr magneton number,^) the 
Brillouin-like magnetization vs temperature curve, '^""'^^^ 
and a large peak of the specific heat at the Curie tem- 
perature [TqY^^^^^ are well explained by a simple local- 
moment model. 

The itinerant vs local-moment behavior in 3(i tran- 
sition metals has been a long-standing problem in the 
metallic magnetism. The Stoner model combined with 
the first-principles band theory yields the Curie temper- 
atures being much higher than the experimental ones; 
6000K for Fe and 3000K for Ni,i6.i7) and does not ex- 
plain the Curie- Weiss law. The theory has been much 
improved by taking into account spin fluctuations. Hub- 
bard'^*-' and Hasegawa^^^ proposed a single-site spin fluc- 
tuation theory (SSF) on the basis of the functional inte- 
gral method. ^''"^^■' The theory interpolates between the 
weak and strong Coulomb interaction limits and ex- 
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plained qualitatively the magnetization vs temperature 
curve as well as the Curie- Weiss law on the basis of the 
band model. The theory however is based on the high- 
temperature approximation, i.e., the static approxima- 
tion which neglects the time dependence of the field vari- 
ables. Therefore the SSF reduces to the Hartree-Fock ap- 
proximation at the ground state, and thus does not take 
into account electron correlations as found in the ground 
state theories. ^^"^^^ 

Kakehashi and Fulde"^^' proposed a variational the- 
ory which takes into account the electron correlations 
at the ground state and reduces to the SSF at high- 
temperatures. They found that the reduction of the 
Curie temperature by a factor of two in the case of Fe. 
Hasegawa'^'^-' developed the same type of theory on the 
basis of the Slave-Boson functional integral approach. 
Although these theories include the ground-state elec- 
tron correlations within the Gutzwiller-type approxima- 
tion, ^^"^^'■^^^ systematic improvement of thermodynam- 
ics is not easy in these approaches. We therefore pro- 
posed a theory called the dynamical coherent potential 
approximation (CPA)'^'^^ which completely takes into ac- 
count the dynamical charge and spin fluctuations within 
the single-site approximations, and clarified numerically 
the basic aspects of the theory with use of a Monte-Carlo 
technique. 

In order to simplify numerical calculations, we pro- 
posed in the next paper, which we refer to 1, more 
analytic theory of the dynamical CPA using the har- 
monic approximation, and verified on the basis of the 
Hubbard model that the dynamical effects reduce the 
Curie temperatures, cause the band narrowing for quasi- 
particle states, and create the '6 eV satellite peak in ex- 
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citation spectra. In our recent paper, which we refer 
to II in the following, we proposed the first-principles 
dynamical CPA. The theory combines the dynamical 
CPA with the Local Dcinsity Approximation (LDA) + 
U scheme^^-* in the tight-binding linear muffin-tin or- 
bital (TB-LMTO)^'') representation. We investigated the 
dynamical effects on the magnetic properties of Fe and 
Ni within the 2nd-order dynamical corrections. Quite re- 
cently, we have improved the first-principles dynamical 
CPA taking into account the higher-order corrections,''^^ 
and clarified the systematic change of excitation spectra 
in 3d transition metal series. 

In this paper, we present numerical results of calcu- 
lations for the magnetic properties of Fe, Co, and Ni 
which are obtained by the 4th-order first-principles dy- 
namical CPA. We investigate the dynamical effects on 
various quantities, and clarify the quantitative aspects 
of the first-principles dynamical CPA on the magnetic 
properties. 

As we have proven in our previous papers,^' ■'^•' the 
dynamical CPA is equivalent to the many-body CPA^°' 
in disordered alloys, the dynamical mean-field theory 
(DMFT)^^"^^) in the metal-insulator transition, and the 
projection operator CPA^^' in excitation problem in 
solids. The first-principles DMFT calculations for Fe 
and Ni at the ground state have been performed by 
Miura and Fujiwara^'"' within the iterative perturbation 
method. The finite-temperature DMFT calculations for 
Fe and Ni have been performed by using the Hamiltonian 
without transverse spin fiuctuations.^^' We present here 
the finite-temperature results for the Hamiltonian with 
transverse spin fiuctuations including the results for Co. 

In the following section, we outline the first-principles 
dynamical CPA starting from the TB-LMTO type of 
Hamiltonian with intraatomic Coulomb interactions, and 
elucidate how to take into account the higher-order dy- 
namical corrections using the asymptotic approximation. 
In §3, we present the numerical results of the densities 
of states, effective potentials, magnetization vs temper- 
ature curves, paramagnetic susceptibilities, and the am- 
plitudes of local moments as a function of temperature. 
We will discuss on the dynamical effects on these phys- 
ical quantities and quantitative aspects of the theory. 
Especially, we will demonstrate that high-temperature 
physical quantities such as the effective Bohr magneton 
number arc qualitatively described by the present the- 
ory. We summarize our results in the last section §4, and 
discuss remaining problems. 

2. First-principles TB-LMTO dynamical CPA 

We adopt in the present paper the TB-LMTO Hamil- 
tonian combined with a LDA-I-U Coulomb interactions 
as foUows.^^) 



H = Hn + H 



(1) 



Ho = ~ ^»^<^ 



^ UlJL' diLa^^jL'a , (2) 
iLjL'a 



iLa 



i m m>m' 

^f^ilmX^ilrn' ^ ^ J^ilm ' ^i/m'j • 

m>m' 



(3) 



Here we assumed a transition metal with an atom per 
unit cell. 6° is an atomic level on site i and orbital L, 
^ is the chemical potential, tiLjL' is a transfer integral 
between orbitals iL and jL'. L = {l,m) denotes s, p, 
and d orbitals. aj^^ {diLa) is the creation (annihilation) 
operator for an electron with orbital L and spin a on site 
i, and hn^rr = (^iLa'^iLa is a charge density operator for 
electrons with orbital L and spin a on site i. 

The Coulomb interaction term Hi consists of the on- 
site interactions between d electrons {I = 2). Uo {Ui) 
and J denote the intra-orbital (inter-orbital) Coulomb 
and exchange interactions, respectively, num (sum) with 
/ = 2 is the charge (spin) density operator for d electrons 
on site i and orbital m. Note that the atomic level e° 
in Hq is not identical with the LDA atomic level e^; 

= €l — dEj^Y)A/'^^iLa- Here UiLa is the charge density 
at the ground state, S^da is a LDA functional to the 
intraatomic Coulomb interactions.'^^' ''^^ 

In the dynamical CPA, we transform in the free energy 
the interacting Hamiltonian Hi into a one-body Hamilto- 
nian with dynamical potential v for time-dependent ran- 
dom charge and exchange fields, using the functional in- 
tegral method. ^^'^^^ Introducing a site-diagonal uniform 
medium, i.e., a coherent potential S into the potential 
part, we expand the correction v — 12 with respect to 
sites. The zeroth term in the expansion is the free en- 
ergy for a uniform medium, .F[E]. The next term is an 
impurity contribution to the free energy. The dynami- 
cal CPA neglects the higher-order terms associated with 
inter-site correlations. The free energy per atom is then 
given by^'*' 



7-CPA = - r^ln J [ny^dC.] e-^^eff 



(4) 



Here f3 is the inverse temperature, = Jy = J± = 

[1-1/(2/+!)] J, X = [/o/(2;+l)-HJi, and ^= (e.,?^;,^.) 
is a static field variable on a site. 

The effective potential £'efr(^) in eq. (4) consists of the 
static contribution Est{^) and the dynamical correction 
term Edyn{$) as follows. 

EM$) = E,,{$) + Edy^{$). (5) 

The former is given as 

mn 

1 ^0 



1 



{Uo-2Ui+j)j2Miy 
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-{2u,-j)ni{if+jie^+jie 



(6) 



Here 5vLa{^) = wl^(O) - SL^(iw„), and = ^2 _^ 
i^Lo-(O) is a static potential given by wl(t(0) = [(?7o ^ 
2J7i + J)n;„(C) + (2i7i - J)n,(4)]/2 - X^^^/S, Sl,(ic^„) 
is the coherent potential in Matsubara frequency repre- 
sentation, and ojn = (2n + l)7r//3. The electron number 
for a given ^ is expressed by means of an impurity 
Green function as 



(7) 



and ni{$^) = X^m^-iC^)- ^^'^ impurity Green function 
Ghcriji^i^n) has to be determined self-consistently. The 
explicit expression will be given later (see eq. (32)). 

The coherent Green function FL„{iuJn) in eq. (6) is 
defined by 



(8) 



Here {Ho)iLajL'a is the one-electron Hamiltonian 
matrix for the noninteracting Hamiltonian Hq, and 

{'S{iu}n))iLajL'a = ^ La ii^^n)5tj S LL' ■ 

The dynamical potential iJdyn(C) in eq. (5) has been 
obtained within the harmonic approximation.^^' ^^'^'^'^^^ 
It is based on an expansion oi E(\y^{^) with respect to the 
frequency mode of the dynamical potential VLaa'ii^v), 
where oj^ = 2^'7r//3. The harmonic approximation is the 
neglect of the mode-mode coupling terms in the expan- 
sion. We have then 



£;dyn(^) = -/ 



Mn 



1 



1) 



(9) 



Here the determinant Di, is a contribution from a dy- 
namical potential VLaa'{ii^iy) with frequency co^, and the 
upper bar denotes a Gaussian average with respect to the 
dynamical charge and exchange field variables, Cm(^w„) 
and ^rna{i(^n) (a = x,y,z). 

The determinant Di/ is expressed as'^^-* 

1^-1 [21+1 

n l[DAk,m) , (10) 

fe=0 Lm=l 



11 

CL^(i/,m) 1 1 

ak^,m) 1 1 

auik(y,m) 1 







1 



a2,^{v,m) 



.(11) 



Note that 1 in the above determinant denotes the 2x2 
unit matrix, a„(j^, m) is a 2 x 2 matrix defined by 

an{iy,m)acr/ = E VLaa"{i'^iy)9La"a"'{i'^n-ii^u) 

a"a"'cr"" 

XVLa"'a""{-ii^u)gLa""a'iii^n) , (12) 



a m' 



ghaa'ii'^n) = [(-FL(iw„) ^ - 5vo) ^] 



(13) 
(14) 



Here CTq. (a = x, y, z) are the Pauli spin matrices. Amm> = 

UoS„,m' + {2Ui-J){l-S„,m'), = Jil-Srnm') (« = 

X,J/), and -B,^,„,, = UaSmm' + Ji^ - Smm')- gLaa'iii^n) 

is the impurity Grecm function in the static approxima- 
tion, {FL{ioJn))aa' = FLa{i^^n)5aa' , and 5vo is defined by 

The determinant D^{k,m) is expanded with respect 
to the dynamical potential as follows. 

D,{k,m) = l + Dl^\k,m) + D^J\k,m) + --- , (15) 

a;i7i---«n7n 

(16) 

Here the subscripts aj and 7i take 4 values 0, x, y, and 
2;, and 



(17) 



Vc{u,m) = -i^B^„„t„/„(iw^)<5;2 , {a = x,y,z) . 



(18) 



Note that the subscript {0:7} of _D|'^^^j(i/, A:, m) in eq. 
(16) denotes a set of (ai7i, • • • , Q:„7„). 

Substituting eq. (15) into eq. (10) and taking the Gaus- 
sian average, we reach 



00 00 



i^dyn(0 = -/3-nn l + ^^i?t' 



■(») 



(19) 



n=l 1^=1 



and 



D 



(n) 



1 



(2/?) 



2;+i 

E E Enn 



Efc^Kfe>^)='^{aj(fe>^)} P m=lfe=0 



.(20) 



Here each element of {l{k,m)} {k = 0, ■ • • , f — 1, m = 
I,-- - ,2Z -|- 1) has a value of zero or positive integer. 
aj{k,m) takes one of 4 cases 0, x, y, and z. j denotes 
the j-th member of the {k,m) block with l{k,m) el- 
ements. P denotes a permutation of a set {{j,k,m)}; 
P{(j, fc,m)} = {(jp, fcp, mp)}. ap-i means a rearrange- 
ment of {aj{k, m)} according to the inverse permutation 
P~^. The coefficient C^„, in eq. (20) is a Coulomb in- 
teraction defined by 



Co 



(a = 0) 

(a = x,y,z) . 



(21) 
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The frequency dependent factors -D|"^^j (i/, fc, m) in eq. 

(20) consist of a linear combination of 2n products of 
the static Green functions. Their first few terms have 
been given in Appendix A of our paper 11.^^^ 

In the calculations of the higher-order dynamical cor- 
rections'^®^ Z)|^'^j(z^, k, m), we note that the coupling con- 
stants = -B^to' = "^(1 ~ Smm') are much smaller 
than Amrn' and B^^, because Uq and Ui ^ J. Thus we 
neglect the transverse potentials, Vx{'^,m) and Vy{v,m). 
The approximation implies that a„(i^, m)^-^ = 0. The 
determinant D^{k, m) in eq. (10) is then written by the 
products of the single-spin components as 



Di,{k,m) = D„^{k,m)D„i{k,m). 



(22) 



Here Dva{k, m) is defined by eq. (11) in which the 2x2 
unit matrices have been replaced by 1 (i.e., 1 x 1 unit 
matrices), and the 2x2 matrices a„(i/, m) have been 
replaced by the 1x1 matrices o„(z/, m)^^. The latter is 
given by 



0,z 

E 

a, 7 



Vail', m)Vj{-l', m)hajaena{l', m), (23) 



(24) 



Here ha^^ = ^a-f + ct(1 — S^-y), and we used a notation 
9L<j{n) = gLaa{i'^n) for simplicity. 

In order to reduce these summations, we make use of 
an asymptotic approximation. 



e„a{y, m) ~ {gLa(n -v)- gLcrinj) , 



(25) 



where q^, = p/i-Kvi. The approximation is justified in the 
high-frequency limit where gLa{n) written as 



9La{n) = - 



+ 



1 



In the asymptotic approximation, we obtain 



• (26) 



xD 





(n-l) 



1 + 



-ghaiiv + k) 



(30) 



Here L>|°^^j^(i^, fc, m) = 1. ni{{a-f}) i 
{ttiTi} pairs such that ai = among the / pairs. When 
there is no orbital degeneracy, eq. (30) reduces to the re- 
sult of the zeroth asymptotic approximation in our paper 

j34) 

In the actual applications we make use of the exact 
form up to a certain order of expansion in D^^\ and 
for higher order terms we adopt an approximate form 
(27). In this way, we can take into account dynamical 
corrections systematically starting from both sides, the 
weak interaction limit and the high-temperature one. 

The coherent potential can be determined by the sta- 
tionary condition SJ^cpa/S'E = 0. This yields the CPA 
equation as^^^ 

{GLa{^,iu„)) = FU^^r.) ■ (31) 

Here ( ) at the l.h.s. (left-hand-side) is a classical average 
taken with respect to the effective potential Eeg{^). The 
impurity Green function is given by 

EE 

Gl<t(^,«Wj) = gLaa{i<^l) + 



is the number of 



5DI 

KL„{iuJl)ST,L„{iLJl) 



■(") 



(32) 



Note that the first term at the r.h.s. (right-hand-side) 
is the impurity Green function in the static approxi- 
mation, which is given by eq. (14). The second term 
is the dynamical corrections, and KL„{i(jJi) = 1 — 

FLa{i'^l)~'^5FLa{ii^l)/5^La{ii^l)- 

Solving the CPA equation (31) self-consistently, we ob- 
tain the effective medium. The electron number on each 
orbital L is then calculated from 



(33) 



{ai + i7i+i---Q:„7„}4, 

Here we wrote the subscript at the r.h.s. explicitly to 

avoid confusion. Note that the values of ai and 7^ are 
limited to or 2: in the present approximation. The spin- 
dependent quantities are given by^^^ 



The chemical potential /i is determined from the condi- 
[y, k,m) . (27) ^[q^ _ X^ii^^)- Here Ue denotes the conduction elec- 



tron number per atom. The magnetic moment is given 

by 



(34) 



D 



(0 

{07} 



^{i^,k,m) = Al}\{aj})f^B^J\u,k,m) , (28) 



In particular, the I = 2 component of magnetic moment 
is expressed as 



A«({"7}) = 



1 (a=t) 
(_l);-"i({«7}) =;) 



(29) 



(35) 



/-I 

B^){v,k,m) = [[{-gLaijy + k) 



3=0 



The above relation implies that the effective potential 
EesiO is ^ potential energy for a local magnetic moment 

. . . 

In the numerical calculations, wc took into account 

the dynamical corrections up to the second order (n < 2) 
exactly, and the higher-order terms up to the fourth order 
within the asymptotic approximation. Summation with 
respect to 1/ in Eqs. (19) and (32) was taken up to = 100 
for n = 1 and 2, and up to = 2 for n = 3, 4. 
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Fig. 1. Up and down d partial densities of states (DOS) of Fc 
in the static approximation (dotted curves) and in the 4th-order 
dynamical CPA (solid curves). The total DOS in the 4th-order 
dynamical CPA is shown by dashed curve. 



When we solve the CPA equation (31), we adopted 
a decoupling approximation to the thermal average of 
impurity Green function^^^ for simplicity, i.e., 



q=± 



1 + q- 



(6) 



^,(ei),»u;„).(36) 



Here wc wrote the static exchange field ^ as (Cz,^!) 
so that the decoupling approximation we made becomes 
clearer. The approximation is correct up to the second 
moment {i.e., (C^)) and arrows us to describe the thermal 
spin fluctuations in a simpler way. 

On the other hand, we adopted a diagonal approxima- 
tions^^ to the coherent Green function at the r.h.s. of eq. 
(31). 



I 



lOJr, 



(37) 



5^L(T(«Wn) - Ae/, 

Here Pi^^ie) is the local density of states for the LDA 
band calculation, and Ae^ = (el — ^'l)5i2- The approxi- 
mation partly takes into account the effects of hybridiza- 
tion between different I blocks in the nonmagnetic state, 
but neglects the effects via spin polarization. 

The CPA equation (31) with use of the decoupling ap- 
proximation (36) yields an approximate solution to the 
full CPA equation. For the calculations of the single- 
particle densities of states, one needs more accurate so- 
lution for the CPA self-consistent equation. For this pur- 
pose, we adopted the following average t-matrix approx- 
imation'''^' '"'S^ (ATA) after we solved eq. (31) with the 
decoupling approximation (36). 

{GLa{S,z,^'i,i0Jn)) - FLgjiuJn) 

Here the coherent potential in the decoupling approx- 
imation is tiscd at the r.h.s., but the full average ( ) 
of the impurity Green function is taken. The ATA is a 
one-shot correction to the full CPA. The coherent poten- 



StJA(iu;„) = Si<.(jw„) + 



tial T,L<j{z) on the real axis z = uj + iS is then calcu- 
lated by using the Pade numerical analytic continuation 
method. Here S is an infinitesimal positive number. 
The densities of states (DOS) as the single-particle exci- 
tations, pLcr{^) are calculated from the relation, 



--Im FLa{z) . 

TT 



(39) 



3. Numerical results 



In the numerical calculations, we adopted the lattice 
constants used by Andersen et. al.,^"^^ and performed the 
LDA calculations with use of the Barth-Hedin exchange- 
correlation potential to make the TB-LMTO Hamilto- 
nian (2). For Fe and Ni, we adopted average Coulomb 
interaction parameters U and the average exchange in- 
teractions J used by Anisimov et. al.,^^'^ and for Co we 
adopted U obtained by Bandyopadhyay et. al.^"^^ and 
J obtained by the Hartree-Fock atomic calculations^^^ ; 
U = 0.169 Ry and 7 = 0.066 Ry for bcc Fc, U = 0.245 
Ry and J = 0.069 Ry for fee Co, and U = 0.221 Ry and 
J = 0.066 Ry for fee Ni. The intra-orbital Coulomb in- 
teraction Ua, inter-orbital Coulomb interaction Ui, and 
the exchange interaction energy parameter J were calcu- 
lated frorn 17 and 7 as Z/q = F -|- 87/5, Ui = U - 27/5, 
and J ~ J, using the relation Uq = Ui + 2 J. 

Figure 1 shows the calculated densities of states (DOS) 
in the ferromagnetic Fe. Thermal spin fluctuations in 
the static approximation broaden the band width for 
each spin band. The DOS with the 4th-order dynam- 
ical charge and spin fluctuations shifts the main peak 
towards the Fermi level and causes a small hump around 
ui = —0.5 Ry which corresponds to the '6 eV satel- 
lite as found in Ni at 6 eV below the Fermi level. ^^"^^^ 
The present calculations with dynamical effects yields 
the ground-state magnetization 2.58 /jb, which is consid- 
erably larger than the experimental value 2.22 /zb- The 
exchange splitting in the present calculations is therefore 
somewhat overestimated. When the exchange splitting is 
reduced, the main peak in the up spin band shifts up and 
the second peak of the down spin band shifts down, so 
that one expects that the main peak around oj = —0.1 
Ry in the total DOS becomes sharper and shifts towards 
the Fermi level. 

The effective potential in the dynamical CPA char- 
acterizes spin fluctuations of the system. Unlike the sin- 
gle band model, '^^^ calculated potential for ferromagnetic 
Fe has the double minimum structure as shown in Fig. 
2. This implies that Fe local magnetic moments show 
large thermal spin fluctuations which change magnetic 
moments in direction. In order to clarify the dynamical 
effects on the effective potential, we plotted the dynam- 
ical contribution -Edyn(^) = £'cff(€) ~ Esi{^) in Fig. 3. 
The dynamical potential of Fe has a saddle point at the 
origin; £^dyn(^) shows the minimum along the z axis, and 
iSc^maximum along the x (y) axis at the origin. It indi- 
cates that the dynamical effects reduce the longitudinal 
amplitude of magnetic moments and enhance the trans- 
verse spin fluctuations. These effects are enhanced with 
increasing temperatures. 

The magnetization vs temperature curves of Fe in vari- 
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Fig. 2. Effective potential for Fe at the temperature T/Tc 
on the — £,z plane. 
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Fig. 4. Magnetization vs temperature curves (M — T), inverse 
susceptibilities (x^^), and the amplitudes of local magnetic mo- 
ments ((m^)^/'^) for Fe in the static approximation (dotted 
curves), the 2nd-order dynamical CPA (dashed curves), and the 
4th-order dynamical CPA (solid curves), respectively. Experimen- 
tal M — T curve^") is shown by + points. 
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Fig. 3. Dynamical contribution to the effective potential for Fe 
at T/Tc = 0.5. 



ous approximations are presented in Fig. 4. We obtained 
the Curie temperature Tq = 2070K in the static ap- 
proximation. It is lower than that of the Hartrec-Fock 
one (12200K) by a factor of 6. The second-order dy- 
namical corrections reduces Tq by 50K as shown in Fig. 
4. The 4th-order dynamical corrections further reduce 
Tc by 90K as compared with the 2nd-order ones, and 
yield Tq = 1930K. The result is comparable to the 
value T = 1900K based on the DMFT without trans- 
verse spin fluctuations,*^^ but still overestimates Tc by 
a factor of 1.8 as compared with the experimental one 
(1040K).^) Calculated inverse susceptibilities follows the 
Curie- Weiss law. We obtained the effective Bohr mag- 
neton number as 3.1 /kb (static approximation), 3.0 hb 
(2nd-order dynamical CPA), and 3.0 /xb (4th-order dy- 
namical CPA), respectively. These values are in good 
agreement with the experimental vahie^' 3.2 /ie- The 
amplitudes of local magnetic moments {m?)^/"^ show a 



Fig. 5. Up and down d partial DOS of Co in the static approxi- 
mation (dotted curves) and in the 4th-order dynamical CPA (solid 
curves). The total DOS in the 4th-order dynamical CPA is shown 
by dashed curve. 



weak temperature dependence and take a value 3.1 /Ub at 
2000K irrespective of details of approximations as shown 
in Fig. 4. It should be noted that the calculated effective 
Bohr magneton number approximately equals to the am- 
plitude of local moment; the Rhodes- Wohlfarth ratio is 
1 in agreement with the experimental fact. 

In the case of Co, the crystal structure changes from 
the hep to the fee with increasing temperature. We 
present here the results for the fee Co. Figure 5 shows an 
example of calculated DOS in the ferromagnetic state. 
The d DOS are split into up and down parts. The d DOS 
in the static approximation are smoothed and are broad- 
ened due to thermal spin fluctuations. The dynamical ef- 
fects reduce the exchange splitting and the band width of 
each spin component. More important difference between 
the static and dynamical cases is that the quasiparticle 
peaks appear at the Fermi level due to dynamical spin 
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Fig. 6. Effective potential of Co at T/Tq = 0.4. 
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Fig. 7. Dynamical contribution to tlie effective potential in Co 
at T/Tc = 0.4. 

and charge fluctuations. Wc also find a smafl hump at 
oj = — 0.55Ry in the up-spin DOS, which is caused by 
dynamical charge fluctuations. 

We present in Fig. 6 the effective potential for fee Co. 
The potential has a single minimum in the ferromagnetic 
state. With increasing temperature the flat part of the 
potential around the origin extends to the negative region 
of and the double minimum structure appears near 
and above Tc- The dynamical corrections to the eff'ective 
potential show a butterfly structure in which the ala in 
the region > is higher than that in < 0, as shown 
in Fig. 7. It indicates that the dynamical eff'ects act to 
reduce the magnetization and to enhance the transverse 
spin fluctuations. 

Calculated magnetization vs temperature curves for 
fee Co are shown in Fig. 8. We find the Curie tempera- 
ture in the static approximation Tc = 3160 K. It is much 
lower than the Hartree-Fock one (12100 K). The 4th- 
order dynamical CPA reduces Tc by 610K, and yields 



3.5 



3 




500 1000 1500 2000 2500 



T(K) 

Fig. 8. Magnetization vs temperature curves {M — T), inverse 
susceptibilities (x^^)i and the amplitudes of local magnetic mo- 
ments ((m^)^/^) for Co in the static approximation (dotted curves) 
and the 4th-order dynamical CPA (solid curves), respectively. Ex- 
perimental data of magnetizations^^^ for fee Co are shown by open 
squares. Dashed line shows a guide for the eye to an experimental 
M — T curve. 

Tc = 2550K. The latter is overestimated by a factor of 
1.8 as compared with the experimental value^^ 1388K. 
The inverse susceptibility follows the Curie- Weiss law 
though it is considerably upward convex. Calculated ef- 
fective Bohr magneton numbers at T/Tq 1.1 are ob- 
tained to be 2.4 /iB in the static approximation and 3.0 
/xb in the 4th-order dynamical CPA. The latter agrees 
well with the experimental value 3.2 /ib. The amplitude 
of local moment hardly change with increasing temper- 
atures, and has a value 2.66 /xb around 2600K. This is 
close to the effective Bohr magneton number 2.55 in the 
local moment model, but is smaller than the experimen- 
tal value 3.2 /iB- 

Calculated DOS for the ferromagnetic Ni is shown in 
Fig. 9. By comparing the DOS with those in the static 
approximation, we find that the dynamical effects reduce 
the exchange splitting by a factor of two and cause a kink 
at the Fermi level corresponding to a quasi-particle state. 
On the other hand, a satellite which is found experimen- 
tally around 6 eV below the Fermi level^^^^'^^ does not 
appear. It indicates that the dynamical corrections up 
to the 4th order with asymptotic approximation are not 
enough to describe the charge fluctuations at low tem- 
peratures T < 500K, though the 6 eV satellite is found at 
ui = — 0.45Ry in the DOS when calculated at high tem- 
peratures in the paramagnetic state (see the thin dashed 
curve and thin solid curve in Fig. 9). 

The effective potential in Ni shows a single minimum 
in both the ferro- and the para- magnetic states as shown 
in Fig. 10, indicating small thermal spin fluctuations 
around the equilibrium points. The dynamical contribu- 
tion to the effective potential in the ferromagnetic Ni 
shows a flat structure with a slope along direction (see 
Fig. 11); -Edyn(^) ~ ^zhes where hes is an effective fleld. 
This implies that the dynamical effects act as an effective 
magnetic fleld which weakens the spin polarization. 

The magnetization vs temperature curves are pre- 
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Fig. 9. Up and down d partial DOS of Ni in tlie static ap- 
proximation (dotted curves) and in the 4th-order dynamical CPA 
(solid curves). The total DOS in the 4th-order dynamical CPA 
is shown by dashed curve. The d DOS in the paramagnetic state 
(T = 2000K) are also shown by thin dashed curve (4th-order dy- 
namical CPA) and thin solid curve (2nd-order dynamical CPA). 



Fig. 12. Magnetization vs temperature curves (M — T), in- 
verse susceptibilities (x~^)i and the amplitude of local magnetic 
moments ((m^)^''^) for Ni in the static approximation (dotted 
curves), the 2nd-order dynamical CPA (dashed curves), and the 
4th-order dynamical CPA (solid curves), respectively. Experimen- 
tal data of magnetization curve^^' are shown by -|-. 
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Fig. 10. Effective potential of Ni at T/Tc = 0.8. 
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Fig. 11. Dynamical contribution to the effective potential in Ni 
at T/Tc = 0.8. 



sented in Fig. 12. We find calculated Curie temperatures: 
1420K (static approximation), 1260K (2nd-order dynam- 
ical CPA), and 620K (4tli-ordcr dynamical CPA), respec- 
tively. These results are much lower than the Hartree- 
Fock value 4940 K. The 4th-order result is in good agree- 
ment with the experimental value^^ 630K. The stability 
of the ferromagnetism in Ni is sensitive to the change of 
the DOS at the Fermi level according to the ground-state 
theory of the ferromagnetism in the low density limit. 
A large reduction of Tq due to 4th-order dynamical cor- 
rections seems to be associated with the instability of 
the ferromagnetism due to the reduction of the DOS at 
the Fermi level. In fact, the band width of the d DOS in 
the 4th-order dynamical CPA is broader than that of the 
the 2nd-order one, and the height of the peak in the 4th 
order DOS is smaller than the 2nd-order one as shown 
in Fig. 9. 

Calculated inverse stisceptibilities in Ni follow the 
Curie- Weiss law in the static approximation and in the 
2nd-order dynamical CPA, but show an upward convex- 
ity in the case of the 4th-ordcr dynamical CPA. The up- 
ward convexity in the high-temperature region is found 
in the experimental data,^^ though it is not the case 
in Co. Effective Bohr magneton numbers calculated at 
T ~ 2000K are 1.2 /xb in the static approximation as 
well as in the 2nd-order dynamical CPA, while 1.6 /kb 
in the 4th-order dynamical CPA. The latter is in good 
agreement with the experimental value^-* 1.6 /^b- Cal- 
culated amplitude of local moment (m^)^/^ slightly in- 
creases with increasing temperature and takes a vahie 
1.97 /iB at lOOOK, which is larger than 1.27 /xb, the value 
in the local moment model. Because of the hybridization 
of the d bands with the sp bands, the d electron number 
in the metallic state (rid = 8.7) is smaller than = 9.4 
expected from a d band model with the strong ferro- 
magnetism, so that the amplitude of the local magnetic 
moment is larger than the value expected from the local 
moment model, yM(0)(M(0) -h 2) where M(0) denotes 
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Fig. 13. Effective Bohr magneton numbers in various 3d transi- 
tion metal alloys as a function of conduction electron number per 
atom.^^^'^^' Large closed circles show theoretical values of Fe, Co, 
and Ni obtained in the present calculations. 



the ground-state magnetization. 

4. Summary and Discussions 

We have investigated the ferromagnetic properties of 
Fe, Co, and Ni on the basis of the first-principles dy- 
namical CPA combined with the LDA + U Hamiltonian 
in the TB-LMTO representation. The dynamical CPA 
takes into account single-site spin and charge fluctua- 
tions. We adopted the harmonic approximation to solve 
the impurity problem in the self-consistent dynamical 
CPA. In the harmonic approximation, we start from the 
high temperature approximation, i.e., the static approxi- 
mation, and take into account individually the dynamical 
contributions from the dynamical potentials with a given 
freqTiency. Wo have performed the numerical calculations 
taking into account the dynamical effects up to the 4th 
order in Coulomb interaction. Calculated effective poten- 
tial shows up the double minimum structure in case of 
Fe, the single minimum structure in case of Ni, and the 
Co is in between the two. 

We found that dynamical effects reduce the exchange 
splitting in the ferromagnetic DOS, suppress the band 
broadenings due to thermal spin fluctuations in the static 
approximation, and cause the '6 eV satellites. In the case 
of Co and Ni, we also found the quasiparticle peaks at 
the Fermi level. 

We calculated the magnetization vs temperature 
curves. Curie temperatures obtained from the 4th-order 
dynamical CPA are 1930K for Fe, 2550K for fee Co, and 
620K for Ni. Although calculated Tq in Ni is close to the 
experimental value 630K, those in Fe and Co are overes- 
timated by a factor of f .8 in comparison with the experi- 
mental values, 1040K (Fe) and 1388K (Co), respectively. 
In the present calculations, the ground state magnetiza- 
tions obtained by an extrapolation are 2.58 (Fe), 1.82 
A*B (Co), and 0.63 iib (Ni), which are considerably larger 
than the experimental values, 2.22 /xb (Fe),''^ 1.74 /Ub 
(Co),^) and 0.62 /xb (Ni).'^) These facts suggest that the 
present calculations tend to underestimate the screening 
of Coulomb interactions, especially at low temperatures. 



It is desirable to take into account dynamical effects more 
accurately in order to suppress the magnetization below 
Tc- One has also to take into account inter-site spin cor- 
relations in order to realize quantitative description of 
the Curie temperature, going beyond the single-site ap- 
proximation. 

We have calculated the paramagnetic susceptibilities. 
Calculated effective Bohr magneton numbers, 3.0 /xb 
(Fe), 3.0 (Ub (Co), and 1.6 hb (Ni), quantitatively explain 
the experimental data as shown in Fig. 13. Experimental 
data of effective Bohr magneton numbers in 3d transi- 
tion metal alloys^^"'^^^ continuously distribute including 
the three points for Fc, Co, and Ni. It is a future prob- 
lem to investigate whether or not the dynamical CPA 
can quantitatively describe these data after extension of 
the theory to alloys. 
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